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Abstract. This paper and the others in the series challenge the standard model of the 
effects of gravitational lensing on observations at large distances. We show that due to the cumu- 
lative effect of lensing, areas corresponding to an observed solid angle can be quite different than 
would be estimated from the corresponding Friedmann-Lemaitre model, even when averaged over 
large angular scales. This paper concentrates on the specific example of spherically symmetric but 
spatially inhomogeneous dust universes, the Lemaitre-Tolman-Bondi models, and shows that radial 
lensing significantly distorts the area distance-redshift and density-redshift relations in these exact 
solutions compared with the standard ones for Friedmann-Lemaitre models. Thus inhomogeneity 
may introduce significant errors into distance estimates based on the standard fl relations, even 
' after all-sky averaging. In addition a useful new gauge choice is presented for these models, solving 

ON . the problem of locating the past null cone exactly. 
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I. INTRODUCTION AND GENERAL ARGUMENTS 



> ' 

-xf : 

Our aim in this series of papers is to show that when analyzing observations at large redshift in the real universe, 
the assumption of an area distance corresponding to that of a best-fit Friedmann-Lemaitre (fl) modelQ — that is, 
one with a matching averaged matter density — may not be a good approximation, even when averaging over large 
angular scales. This claim is important because of the ubiquitous use of fl models in studies of number counts versus 
redshift and area distance versus redshift. 



The first paper in the series (paper I) gave general arguments for this thesis. It was explained there that 



whilst the Dyer-Roeder distance |3j is generally regarded as a good approximation for ray bundles moving between 
?H . high-density clusters of matter, resulting in a de-focussing relative to the comparable fl model, matter moving near 
or through higher density regions is more focused than in the fl model, so resulting in a compensating effect. It is 
J> \ commonly believed |4|,|{| that an exact cancellation between the two effects takes place when one averages over large 
angular scales containing both high and low density regions — so the fl area distance is the correct one on these 
scales. However, it was pointed out in Q that, after caustics have formed through the focussing of light rays by the 
high density regions, these light rays too are rapidly diverging so that sufficiently far down the past light cone all light 
rays will be diverging relative to the comparable fl model. Consequently there is good reason to question the general 
opinion in this matter — we find that in fact 'shrinking' takes place.0 It is often argued that photon conservation 
excludes this possibility, but these arguments are based on incorrect assumptions about the geometry. 



*email: nazeem@maths.uct.ac.za 
^email: bruce@stardust.sissa.it 
"'"email: cwh@maths.uct.ac.za 
^email: ellis@maths.uct.ac.za 

1 We follow recent moves to use 'Friedmann-Lemaitre' to describe the dynamics of and the application of the Einstein field 
equations to the standard model, and 'Robertson- Walker' to describe its metric and geometry. 

2 The term shrinking |j] describes the following. If at some redshift the area distance is greater in the inhomogeneous universe 
than in the best-fit or background fl model, scales at that redshift will be underestimated if the fl area distance is used to 
calculate the size of objects from their measured angular size, and so will have 'shrunk' relative to their actual scales in the 
more realistic model. Put another way, if a sphere at given z about the observer has a larger total area in the inhomogeneous 
model, an object of fixed size will subtend a smaller solid angle on the sky, and so be shrunk. (See [[jj for more discussion.) 
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Given this feature, even when there is only weak lensing|] and caustics have not formed, it is not obvious that the 
FL area distance-redshift relation is an accurate description. 

The aim of this paper is to construct an exact inhomogeneous model and its fl approximation, and compare the 
area distance-redshift and density-redshift relations in the two. To do so we examine Lemaitre-Tolman-Bondi (ltb) 
spherically symmetric dust solutions where exact integrations of the field equations are available for the past 

light cone of observers at the central position. Although the lensing that occurs for this central observer is purely 
radial, we find the inhomogeneity has a tangible effect on observational relations .[] In the particular case examined, 
because of the spherical symmetry about the observer, the effect in any specific direction will not be compensated 
by an opposite effect in another direction — on the contrary, the effect is uniform because it is the same in all 
directions, and will occur on large as well as small angular scales. However unlike the previous paper Q|, this will not 
be associated with the formation of caustics. 

The failure of the FL-like area distance-redshift assumption will thus have been shown to occur for observers at 
a very special position in a family of high-symmetry space-times. Undoubtedly this is not the kind of situation the 
authors of the papers mentioned above, claiming the effect does not occur, had in mind. However the usual statements 
of this effect contain no clauses that exclude this situation: the result is supposed to hold in all cases when the 'lumpy' 
universe is reasonably close to an fl model, and observations are averaged over the sky; their arguments do not exclude 
the situation envisaged here where the observer is at the centre of a spherically symmetric inhomogeneity. Indeed, a 
spherically symmetric model may be regarded as describing data that has been averaged over the whole sky, but not 
over distance. 

Our example thus confirms the claims of paper I, in the setting of particular exact inhomogeneous solutions of 
the Einstein Field Equations. It does not generically establish the magnitude of the effect, precisely because the 
high-symmetry geometry considered here precludes formation of caustics and the consequent fractal-like structure of 
the real light cone. Paper III [|| will provide confirmation of the overall shrinking effect due to caustics, and attempt 
more realistic estimates of its magnitude than those given in paper I, which used simple analytic formulae for this 
purpose. 

In developing the results of this paper, we solve one of the problems that has made analysis of observations in 
Lemaitre-Tolman-Bondi solutions difficult, namely the problem of precisely locating the past light cone of the chosen 
central event P, by use of a special choice of radial coordinate that ensures a very simple form for the past light 
cone of P in these inhomogeneous space-times. This technical development has other uses in terms of analysing 
observational relations in these models.^] Our approach is complementary to that of Kurki-Suonio and Liang |p^| , who 
did numerical calculations of observational relations in a hyperbolic LTB model derived from an Q = 0.1 fl model 
plus some overdensities. 

II. PROGRAMME 

We select the simplest inhomogeneous fluid solution of the Einstein Field Equations; the Lemaitre-Tolman-Bondi 
(ltb) model which is spherically symmetric, but radially inhomogeneous, with a dust equation of state. 

The question we are raising is whether the area of an averaged wavefront we receive at our observatory in an 
inhomogeneous universe is the same as the area of a wavefront in a smoothed version of that universe. To clarify this 
issue our strategy is to: 

• Select the most natural generalisation of the Einstein- de Sitter models commonly used in studies of observations 
(i.e. a parabolic ltb model) and describe data on the null cone. 

• Find the fl limit of this inhomogeneous universe in an appropriate coordinate system. 

• Average the lumpy universe in a natural way and fit it correctly to a fl model. 

• Compare area distances in the lumpy universe and its smoothed average. 



3 We refer to 'lensing' when the light rays are different than they would be in the background fl model. Thus weak lensing 
implies a change of apparent positions and alteration in the usual distance relations, usually combined with image distortions, 
but not necessarily multiple imaging. 

4 Radial lensing is a spherically symmetric distortion of the null cone compared with an fl model, resulting in a uniform delay 
of the wavefront. There is no image distortion, no dependence of magnification or time delay on direction, and no multiple 
imaging, but our results show its effects are observable. 

5 For a slightly different analysis of ltb spacetimes, based on null cone coordinates, see flipf , and for a consideration of 
observations away from the centre of symmetry see JL1|. 
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A. The Inhomogeneous Model 



1. The Integrated Field Equations 

We choose the parabolic ltb model which is the natural generalisation of the f2 = 1 dust fl model. This model is 
characterised by the mass within comoving radius r, M(r), and the 'bang-time' function ^(r) describing the locus 
of the initial spatial hypersurface (that is, the local time of the big bang). The two arbitrary functions represent a 
physical freedom and a coordinate freedom, e.g. t B (M) and M(r), respectively. 

In normalised comoving coordinates the metric after solving the off-diagonal EFE is 

ds 2 = -dt 2 + {R'fdr 2 + R 2 (d9 2 + sin 2 9d<S> 2 ) (1) 

where R'(t,r) — dR(t,r)/dr. 

The time curves are irrotational and, for comoving dust (p = 0), are necessarily geodesies because of momentum 
conservation. The spatial sections are flat because if we choose r = R(t ,r) then R'(t ,r) = 1 and we find that the 
3-spaces have metric da 2 = dr 2 + r 2 dfl 2 and hence are flat. 

The areal radius, R = R(t, r) in the Lemaitre-Tolman-Bondi metric, is the area of the intersection of our past null 
cone with past spacelike time surfaces (in this case spheres). In the parabolic case R is given explicitly by the solution 
to the equation of motion 



obtained from the 11, 22 and 33 components of the EFE, where ' denotes the derivative with respect to t; i.e. 



R(t,r) 



9M(r) 



(* -**(»■)) 



1/3 



(3) 
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where t is cosmic time whilst M and t B are both functions of coordinate radius r only. It follows immediately that 
R'(t,r) = [M'(r)(t-t B (r)) 2 -2M(r)(t-t B (r))t' B (r)] (4) 

where ' denotes the derivative with respect to radial coordinate r. 
The 00 field equation gives 

/ \ M'(r) 

i7Tp{t > r) = RHt,r)R'(t,r)- (5) 



2. The Solution on the Null Cone 



Since we are interested in observations on the null cone we must project onto it by specifying the unique relationship 
between r and t. On radial null geodesies, ds 2 = = d6 2 = d$ 2 , so from ([!]), if the past light cone of the event 
(t =to,r = 0) is given by t = t(r), then that light cone is described by 

dt = ±R'(i(r),r)dr . (6) 

The coordinate freedom in the ltb metric is a rescaling of the radial coordinate r — ► f = f(r). If we choose r so that 

R'(t(r),r) = l, (7) 

then on the past light cone dt — —dr, so that the incoming light rays at the event (t = to, r = 0) are given by 

t(r) =to-r. (8) 

So this gauge choice, in contrast to other work done on observations in the LTB model, locates the null cone of the 
observer at one instant of time, in its simplest possible form, and makes our programme analytically solvable. On 
this light cone, putting (g) in (|j) and using (Q), 
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R'(t(r),r) = 



2R 2 {t(r),r) 



[M'(r)T(r) 2 + 2Af (r)r(r)r'(r) + 2M{r)r{r)] = 1 



where we have defined 



r(r) = t - r - t B (r) . 



(9) 



(10) 



The function r(r) can be interpreted as proper time from the bang surface to our past null cone along the particle 
worldlines. We can set to to be the time since the bang at the observer (r = 0) by choosing is(0) = (so t(0) = to). 

It is important to realise that evaluating R'(t,r) on the null cone t = t(r) is not the same as differentiating 
R(r) — R(t(r), r) with respect to r. In fact, by evaluating (||) on the null cone, R is given by 



R = R(t(r),r) 



9M(r) 



r(r) 



1/3 



which means that its derivative is given by 
dR d 



\R(t(r),r)} = , 

dr dr 1 (K h n 2R 2 (t(r),r) 



[M'(r)r(r) 2 + 2M(r)r(r)r'(r)] 



Combining the above equation with the constraint @ gives a first order differential equation for R. 

dR{r) 3M(r)r{r) 



dr 



= 1. 



R(r) 



(11) 



(12) 



(13) 



In summary, with our choice of coordinates we have recast the flat LTB model in a form that allows us to locate the 
past null cone with ease. This has left us with one physical freedom to choose an arbitrary function of r. We could 
choose r (or M) and substitute ( fTl| ) into (13). Solution of this differential equation would determine R and thus any 
other quantity. If we instead decide to choose R, that is, the area of the wavefront, then the model is trivially and 
fully specified by (O) and @. It follows that 



and 



M{r) 



-(r) 



R(r) 



2R(r) 



/ I 



dR{r) 
dr 



dR{r) 
dr 



(M(0) = 0) . 



(14) 



(15) 



To obtain the results of the next section, we will fix the geometry by choosing R. This will then determine r(r) 
(equivalently ^s(r)) and M(r) by the above two equations. The flat LTB model will thus be fully specified in these 
coordinates and one could then propagate the data off the null cone by the comoving assumption. 
The density on the null cone p(r) is found by evaluating (S) on the null cone: 



47r/5(r) = 



M'(r) 
R(rf 



and its value at the origin depends on the time as characterised by the Hubble constant ;[] 

~ 4?rp Hi 2 

^l^ 1 => "=4^' t0 =3lfo- 



(16) 



(17) 



Since measurements of the Hubble constant are taken at about z < 1, we can take this to determine the age of the universe, 
to, at the central observer. 
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3. Redshifts 



It is of some importance that we state the relevant quantities in terms of redshifts. To do this, we use the fact 
that in the geometric optics limit, for two light rays emitted on the worldline at r em with time interval St ern = 
t + (r em ) — t~(r em ) and observed on the central worldline with time interval 5t b = t + (0) — * _ (0) 



l + z = 



St 



oh 



St,. 



(18) 



The past radial null geodesies are given by 

dt = -R'(t,r)dr, 

so for an observer on a nearby worldline, the time interval changes by 
d(St) = dt + -dr = [R'{t~,r) - R'(t + ,r)]dr 

Thus 

din** = [R'(t,r)]dr 



-[R'(t,r)]5tdr. 



which means that the rcdshift, given by (18), is 

ln(l + z) 



R'(i,r)dr 



where t(r) is the equation of the null cone Q To calculate R'(t,r) we differentiate (||) with respect to t 



it' 



R 



t' B , M' 



t-t B M 



R = 



AM 



-l 1/3 



3(*-* B ) 



1 


'AM' 


1/3 


' M' 1 + r' 


3 


3t 




M T 



Since i{r) = to — r when we choose R' = 1 on the null cone, R'(t, r) is given by 

After some manipulation of the above expression substituted into ( |19|) , we find that 

_ /4M\ 1/3 1 r- /4M^ 1/3 



ln(l + z) = 



V 3r 



3 Jo V 37 



dr . 



(19) 



(20) 



(21) 



(22) 



Using (H) and (J15|) this equation may be written as 



dR 



dR\ , . , 



(23) 



so we can now determine the redshift-area distance relation. 



7 The standard formula 1 + z = {v^k^em/iu^k^ob is not useful in this gauge since — (R' , —1, 0, 0) = (1, —1, 0, 0) is not 
geodesic, though it is tangent to the past null cone. 
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B. The Friedmann-Lemaitre Limit 



The characterisation of the FL limit is that the bang time surface is simultaneous. So ts(r) = ts E 
from whence 



^M FL (r)(t-i BFL ) 2 



1/3 



2RI 



constant; 



(24) 



The freedom left here in M FL (r) is just essentially the coordinate freedom, corresponding to the freedom of choice of 
r. The above relations determine the fl density 



Pfl(*) = 



1 



6ir(t - iB FL ) 2 



(25) 

FL = 0. We do not have a 



which is spatially homogeneous as required, unaffected by M FL (r). It is usual to set ts 
freedom to rescale the density by a constant because this is the critical density case. 

As we would eventually like to compare our LTB model as chosen above to an underlying fl model, it is appropriate 
to write the fl limit in the same kind of coordinate system. Consider light rays coming in to the event (t = ti, r = 0) 
in a fl model. When we choose coordinates for which R' Flj (t, r) = 1 on the null cone, the past null cone can be located 
by t = t\ — r — ts FL =t\ — r. (We use t\ rather than to here, as we will need to distinguish ltb and fl values later 
on.) As a limit of the flat ltb model in these coordinates, the fl form of M(r) is obtained from setting r = t\ — r in 
@. This yields 



ti /S - (ti - r)V3 

.1/3 



M FL = 6 
i? FL = 3 

We note that this in conjunction with (|24|) implies that 



(M FL (0) = 0) 



if-(ti-r) 1 / 3 (ii-rf/ 3 (i? FL (0) = 0). 



1/3 



(h-r) 



1/3 



.2/3 



RL. = 



3M' 



2i? FL 



FL t 2 = 



f 2/3 



The RW metric that results is, from ( P8[) and (|l|), 

1 



ds z 



-dt 2 + i 4 / 3 



■dr 2 + 9 



t{ /3 -(h-r) 



(il - r) 2 / 3 



(26) 
(27) 

(28) 
(29) 



These coordinates are singular at the particle horizon, r = t\ (when the past null cone of t = ti runs into the initial 
singularity). Thus they are valid for < r < t±. The fl redshift-distance formula can be obtained by inserting the 
fl forms of M(r) and r{r) into equation (^2|). That is 



z(r) 



2/3 



- 1 



-CO = ti 



(1 + zfl 2 - 1 



(l + z)3/ 2 



(30) 



C. Averaging and Fitting 

We want to compare and contrast total areas of wavefronts at given redshifts of an inhomogeneous model to that 
of the corresponding fl model of density equal to the inhomogeneous density perfectly smoothed. This must be done 
with respect to the inhomogeneous metric because physically the smoothing does not occur. 

Perhaps the crucial part of our analysis is that we ensure that we compare with the fl model with the correct 
average density. We define the average or background fl model to be the one that matches on at the particle horizon 
where r — 0, r = rs, using the Darmois-Israel boundary conditions |l3|,[l4j]. Matching first and second fundamental 
forms of this timelike (comoving) boundary surface S gives 

-^ltbIx; — Rfl\-£ 

and the background model must be parabolic if the inhomogeneous one is; or vice versa. 



G 



The matching must hold over all of S; that is, at all times — so 



(32) 



and thus, by 



Mltb| s = M FL | E 



LTBlj] 



(33) 







Thus it is sufficient to match the masses at £, and synchronise the starting times (bang times) when R 
i? FL | s . In general, we do not expect the fl radial coordinate on £ (r FL |s) to be the same as the LTB one there 
(^ltbIs — r s) since the coordinate condition R' = 1 holds on the null cone, whose locus is model dependent. 

For a parabolic ltb model with metric ([l]) and density given by (]|), the background density p FL is the same as that 
obtained by integrating over constant time slices. 



2-rr />7r />r£ 



JO JO 



2tt />7r />rx: 



/•*..., = (/ / / fjH drdSd*] I (J J ' [~ y/>jdrdM*\ 



3 M(r E ) 



4tt [i?(r s ,t )] 3 
1 

6Tr(t -t B (rvW 



(34) 



where equation (g) was used. 

One important point that must be made here is that a covariant averaging procedure does not exist as yet. We 
have used here an averaging method which is 'natural' for the comoving synchronous coordinates which lead to a 
3 + 1 foliation of spacetime. However, the same model in different (for example observational) coordinates would 
suggest a different averaging procedure which could conceivably yield different results. Therefore the claim Q that 
the wavefront areas obtained in the inhomogeneous model and the averaged model are the same already seems highly 
unlikely. 



III. RESULTS 



We use geometric units such that G = c = 1. If we choose a unit of time Tq seconds to be 1 geometric time unit 
(gtu), then the geometric units of length, mass, density, etc. are fixed by 1 glu — Lq = cTq metres, 1 gmu = Mq = 
(c 3 /G)Tc kg, 1 gmu glu~ 3 = pq = (1/G)Tq 2 kg m -3 . For the purposes of this paper, we want units suitable to 
cosmological scales, so we specify that one cosmological time unit, 1 ctu, is ten billion years — of the order of the age 
of the universe. This gives us 

Cosmological Geometric Units 





Time 


Length 


Mass 


Density 


Cosmological 


1 ctu 


1 clu 


1 emu 


1 emu clu~ 6 


SI 


3.156 x 10 1Y s 


9.461 x 10 2b m 


1.275 x 10 ba kg 


1.505 x 10~ 2b kg/m A 


Astronomical 


10 Gyr 


3.066 Gpc 


6.409 x 10 22 M Q 


1.505 x 10~ 28 g/cc 



The first subsection ( IIIB ) gives a very simple model which satisfies the criteria for a reasonable cosmological model 
(with the classical Copernican principle dropped) and which provides a proof that there exist physically reasonable 
density behaviours which lead to a nonzero magnification or shrinking. It is obvious that averaging over the sky will 
not remove this effect since the model is already spherically symmetric. The second model (section IIIC) does the 
same, but is smoother at the origin and displays interesting behaviour in redshift space, and the third model (section 
HID) is exactly homogeneous beyond a given comoving radius. 
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A. Form of Perturbation and General Results 



It is not easy to choose a form of area distance function for the inhomogeneous model which results in reasonable 
physical behaviour. So instead we choose it in the form of a 'perturbation' from a flat Friedmann model; that is, 



R(r)=R^(r)(l + 8(r)) 



(35) 



where, from fl27|), i?" L (r) = 3\t u 1 / 3 — (t u — r) 1 ^ 3 } (t u — r) 2 ^ 3 is the area function of an underlying fl model of age 
ti = t u . (This 'underlying' fl model is a mathematical device with no physical significance. It can not be considered 
a background or average model since we have not restricted 8(r) to average out to zero in any sense.) In principle, 
one should choose a density function and then determine the area distance function from it or risk the possibility of 
assuming the result. However, if we can show that the above choice of R leads to a density profile with reasonable 
physical behaviour, this would suffice — since if we had initially chosen that density function, it would lead to an R 
as chosen above. We will show that this is indeed the case and also indicate that the model is free of shell crossings. 8 
Obviously, for 8(r) smooth and finite, R(r) is zero at the same places as i?" L (r), i.e. at r = and at r = t u . For 
this form of perturbation, in terms of the convenient parametrisation v = r/t u , we find 



X 



M = 



1 



(l-v) 1 / 3 



1 



-t u X(l-v)(l 



S)[2X(l + 8)-S- 
2t u X(l-v)(l + 8) 



3t u X(l - v)Sf 



t B 

dln(l + z) 

dz 



[2X(1 + 6)-5- 3t u X(l - v)8' 
to-r-r 

2X(l + 8)-8-3t u X(l-v)S' 



9[t u X(l-v)(l + S)}^ 



{2X(3 + 48) (I + 8)- 8(1 + 8)- 3t u X(l - v)(5 + 68)8' 



+ 36t u X 2 (l - v)(l + 8)8' - 9t u 2 X 2 (l - v) 2 [2(l + 8)8" + 8' 2 ]} 
= {4X(1 + 28) (1 + 8)- 8 2 - 6t u X(l - v)(2 + 38)8' 
+ 36t u X 2 (l - v)(l + 8)8' - 9t u 2 X 2 (l - v) 2 [2(l + 8)8" + 8' 2 }} / [m u X{l -«)(! + 8)] . 



(36) 
(37) 

(38) 
(39) 

(40) 
(41) 



If (5(0) ^ we find the unphysical limits r(0) = and p(0) = oo. Thus we set 8(0) — 0, obtaining the following 
limiting values: 



M(0) 

r(0) 

87rp(0) 
d\n(l + z) 



2(1 - 3t u S'(r)) 2 r 3 



= 



t,. 



dr 



and 



M(t u ) 
r(t u ) 

8np(t u ) 



1 - 3t u 8'(0) 
4(1 - 3t u 8'(0)) 2 
3t u 

2(l-3t u 8'(0)) 
3t u 

6t u (l + 8(t u )) 3 

- (*u-r)U t „ = o 

4(3 + 48(r)) 
' 9(t u - r) 2 



oo 



(42) 
(43) 
(44) 
(45) 



(46) 
(47) 

(48) 



8 The necessary and sufficient conditions for there to be no shell crossings anywhere or at any time in the evolution of a flat 
model with R' > are that M(r) be an increasing and ts(r) a decreasing function. They were found (for all ltb spacetimes) 
by Hellaby and Lake Jl5| . 
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dln(l + z) 



dr 



2(l + 2<J(r)) 



3(i» - r) 



(49) 



(The limits for the background fl model are obtained by setting S = = 8' = 5" and replacing t u by if,.) From 
numerical experimentation we concluded, in order to avoid shell crossings, that S(r) must remain sufficiently far away 
from zero over most if not the entire range of r > 0, and certainly near r — t u . We want the proper time from the 
bang surface to the null surface on the central worldine to be the 'true' age of the universe; that is, we want it to be to, 
the time at the origin of the ltb model. By setting r(0) = to in (ff3[), the age of the underlying model is determined 



*u = to(l-3ti 6'{0 j) 



(50) 



The parameter t u is the r-coordinate value at which the null cone of the LTB model intersects the bang. We will 
average quantities on this scale; that is to say, we shall take = t u . We match this inhomogeneous universe to a 
flat fl model at the surface rs by equating the masses and bang times at that point. This then determines the time 
ii = if, in the background fl model which we will use for our comparison. At r — t u , R — and ( f46| ) shows that at 
this point, 

M(t u ) = M« L (i„) (1 + S(t u )) 3 = 6t u (1 + S(t u )) 3 . 

In the background fl model the value of the mass at £ is 6if, and this is what we have to match the inhomogeneous 
mass to. This gives us a value of the age for the background flat fl model of 



t b = t u (1 + 5(t u )Y 



(51) 



B. A Regular Model which Exhibits Shrinking and Magnification 

The following simple example is physically well behaved, being free of shell crossings at all times in its evolution 
for r <t u . Since t' B =fi at the origin, the model is not as smooth there as one would like, but there are no physical 
problems. We choose 8(r) for our first model, LTB1, to be 

o{r) = — sin . (52) 



5 V i 

When we set r(0) = to = 1 then 



25 

and the age for the background model FL1, after matching the masses, is 

sin0.87r x " 



12tt\ 

t u = to 1 + — (53) 



*6=*u^l —J ■ (54) 

The calculation of the redshift was done by a numerical quadrature of (^3|). Figs. [I] and ^ compare the dependence 
of area distance and density on redshift, R(z) & p(z), in LTB1 with the corresponding functions in the background 
model FL1. 
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FIG. 1. A plot of area distance against redshift R(z) & Rfl(z) on the past null cone of the inhomogeneous model LTB1 and 
the corresponding background model FL1. The units of R are cosmological length units, and all the figures use base 10 logs. 
This shows that there are systematic shrinking (R > Rfl) and magnification (R < Rfl) effects due to purely radial lensing, 
which obviously cannot be removed by averaging over large angular scales or even the whole sky. Effects of true lensing in a 
more realistic universe would be imposed on top of this. 



<Q- 




FIG. 2. The density of matter p(z) & pfl(z) against redshift on the past null cone in the models of FIG. |l[ LTB1 and its 
corresponding background model FL1. The units are cosmological density units (emu clu~ 3 ). When comparing with FIG. [j], 
we see that roughly speaking, magnification occurs for obects in or just beyond an overdense region, and shrinking occurs for 
objects in or just beyond an underdensity. 
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It is important to plot these quantities in terms of the observable quantity z for two reasons. First of all, in the 
transformation r — > z, the possibility exists that the area distances of the flat and inhomogeneous models might 
transform into each other. Secondly, under certain circumstances the redshift becomes disordered with distance and 
unexpected behaviour might occur, as the following model illustrates. 



C. A Regular Model with Multivalued Observations 



As an illustration of how different the physical quantities plotted against radial coordinate r as opposed to those 
same quantities plotted against redshift z may appear, we present here an ltb model for which the redshift becomes 
disordered with distance at some points and then ordered again at later points. 

The universe is chosen as above but with a 'perturbation function' of 



S(r) 



1 . /0.757rr 

- - sin 

4 V U 



1 
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7T 



1 — cos 



(55) 



This model, which we call LTB2, is also free of shell crossings at any time for r < t u and has a completely smooth 
and regular origin (where t' B = 0). Setting r(0) = to = 1 once again gives 



tu — to 

and the appropriate background model FL2 has age 



9tt 
16 



(56) 



h = t u I 1 — 



sin 0.757T 



(57) 



which arc different from those of LTB1. Figs. EM show R and p against r and against z for LTB 2 and FL2. 




FIG. 3. The area distance functions R(r) & Rfl{t) on the past null cones of the second model, LTB2, and its corresponding 
background model, FL2, given in cosmological units. The horizontal variable is v = r/t u or r/tt for LTB2 and FL2 respectively. 
The physical behaviour of LTB2 as plotted against v appears normal. 
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FIG. 4. The densities on the past null cone for LTB2 and FL2, p(r) and Pfl(t), in cosmological units. Again, LTB2's 
inhomogeneous profile vs v appears quite acceptable. 




log (1+z) 



FIG. 5. Area distance against redshift R(z) & Rfl{z) for LTB2 and FL2. The interesting point to note is that at some 
redshifts the area distance in LTB2 is multivalued. For log 1() (l + z) > 1 the graph looks very much like FIG. [l]. 



12 



4 - 



3 - 



2 - 

<CL 

o 1 - 

- 



-1 - 



-2 

0.2 0.4 0.6 0.8 1 

log (1+z) 

FIG. 6. The densities, p(z) and pfl(z), against redshift for LTB2 and FL2. Note the quaint 'looping' behaviour. For 
log 10 (l + z) > 1 the graph looks very much like FIG. ^| 

This model provides a good illustration of why one has to be careful in ascribing physical behaviour in a certain 
coordinate system. Viewed as functions of r, R and p have fairly standard behaviour, but viewed in terms of the 
observable quantity z, the density and area distance become multivalued. Hence, three objects with the same intrinsic 
luminosity located at different distances appear at the same z, with three different apparent luminosities (or area 
distances) . 

Our numerical experiments indicate that the redshift on the light cone is most sensitive to perturbations in the 
vicinity of the maximum in R(z). All our models in which dz / dr became negative did so in this region. The looping 
behaviour in the R vs log(l + z) plot occurs when the maximum and minimum in the log(l + z) vs r graph bracket 
the maximum in the R vs r graph. Similarly, perturbations more easily generate a maximum and minimum in the 
p(z) near the maximum in i?,, hence the loop in that graph. 

D. A Swiss-Cheese Type Model 

In order to emphasise the effect of inhomogeneity on the past null cone, we consider a model which is homogeneous 
on the large scale. We construct a local inhomogeneity which is properly matched onto an exact fl model at some 
comoving radius, r — rj, somewhere between us and the particle horizon. In other words, the outer regions are exactly 
homogeneous. This is a simple model of a universe which is FL on some averaging scale less than the horizon size, 
and is seen by an observer located in an overdensity. 

In particular, for LTB3, we set rj = 0.85t u = 2.352, we retain the inhomogeneous LTB2 model for the interior, 
r <rj, and for the exterior, r > rj, we use the matching fl model, FLJ, with the same mass and bang time as LTB2 
at the boundary Obviously the background model FL3 that we fit at the particle horizon is identically the one 



9 In order to make the radial coordinate continuous at this matching, we use Rflj = 3[t^ 3 — (t2 — r) 1 ^ 3 ](t2 — r) 2 ^ 3 , Mflj = 
6[4 /3 - (t 2 - r) 1/3 ] 3 , where t 2 = T LT B(r 3 ) + r 3 and tlj 3 = T^B^Ti) + (^LTi3(r J )/6) 1/3 . These differ from the purely fl 
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matched at fj, FLJ, but written in the form (|2(])-(^9]) — i.e. with a different r coordinate. 




log (1+z) 



FIG. 7. Area distance against redshift R(z) & Rfl(z) for the swiss-cheese model LTB3 and its corresponding background 
model FL3. The boundary is located at log 10 (l + Zj) = 0.2635, with the model being identical to LTB2 for smaller z, and 
exactly fl for larger z. Though the background and matched fl models are identical, the R(z) curves do not match up in 
the exterior region — there is a constant displacement parallel to the log(l + z) axis. This emphasises the fact that the z 
integral is path dependent, and so the value of z at any given area distance is affected by the warping of the null cone in the 
inhomogeneous interior. 



expressions because our definition of r is also path dependent. It turns out that tb =tz- 
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FIG. 8. Density against redshift p(z) & pfl{z) for the swiss-cheese model LTB3 and its background model FL3. The density 
jumps at the boundary, log 10 (l + 2j) = 0.2635, but this is not unrealistic. (Within the constraints of a parabolic model, it could 
only have been made continuous at the expense of creating shell crossings.) As with the area distance, the two p(z) curves do 
not co-incide in the fl exterior, because the shape of the null cone is determined by the spacetime geometry it encounters along 
its path. The boundary could have been put at any z value, and it is evident from FIGS. |E| and ^ especially that in general the 
standard fl relations are not recovered in the fl exterior. 

Figs and ^ show clearly that the null cone relations R(z) and p{z) in this model (LTB3) do not return to the 
standard fl ones in the fl region beyond the boundary, z > Zj . This is because the redshift z is path dependent and 
the inhomogeneous interior distorts the null cone. After passing through the inhomogeneity, the null cone, and thus 
the value of log 10 (l + z), are displaced relative to a model which is fl throughout. This holds generally for any choice 
of inhomogeneous interior and boundary position. 



IV. CONCLUSIONS 



The general belief that photon conservation implies that the total area of an incoming wavefront must be the same 
as in the background, matter-averaged, fl model is contradicted by our results. The functions R(z) and p(z) in an 
inhomogeneous model differ from their standard fl forms, because the redshift z is path dependent and the null cone 
is warped. Thus, even if the inhomogeneous model is exactly homogeneous beyond some radius, the standard FL 
forms are not recovered. This means that a universe which is fl on average will not in general present standard FL 
forms of R(z) and p(z) to an observer. 

The spherically symmetric model used here is simple but effective, since averaging over direction cannot change the 
results. In more realistic models of the lumpy universe this effect will still be present, and we expect full gravitational 
lensing to occur, resulting in more significant deviations from the FL formula. 

This investigation used a parabolic LTB model, where the areal radius R is also the area distance of the 2-sphere 
wavefronts of the past null cone. The density in the ltb model is averaged to give a background Einstein-de Sitter 
(il = 1) model, and it is tested against this model. Although there exists no covariant way to perform this averaging, 
we use the 'natural' one defined by the use of junction conditions, here equivalent to the one used in astrophysical 
problems: that is, averaging by integrating the rest mass and proper volume on constant time slices. More importantly, 
the functions R(z) and p(z) in the examples are obviously perturbed away from the standard fl ones, so that no one 
FL model can give the same R and p values as an inhomogeneous one at a variety of z values. 

The results show that it is quite easy to have areas in the inhomogeneous models which differ significantly from 
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areas in the background, matter-averaged fl model. The result may either be shrinking (the background fl area 
distance underestimates the real area distance at that redshift) or magnification (background fl area distance is an 
overestimate). The presence of loops in the R-z and p-z graphs is analogous to the well known 'finger of God' effect 
familiar in redshift maps of the galaxy distribution. 

Our results and conclusions generally agree with those of Kurki-Suonio and Liang [jl2| , who calculated obervational 
relations numerically in 4 hyperbolic LTB models out to z = 0.5. They generated mild and strong deviations from the 
FL observational relations, with 'bang time' inhomogeneities and 'areal radius' inhomogeneities having opposite effects 
on the oberved matter distribution in redshift space, compared with the present day distribution. Their method did 
not permit the redshift to be disordered with distance. 

Whilst the major aim of this paper has been the above thesis, the choice of radial coordinate which locates the 
observer's null cone will be of use in future analyses of observations in these isotropic dust models. 

An important caveat is that since the LTB model does not allow for formation of caustics in the null cone of 
the central observer, it cannot be considered a useful model for obtaining quantitative 'real world' results. Rather 
this paper should be viewed as a proof that even purely radial lensing distorts the area distance-redshift relation 
significantly. If the observer moves away from the central position, then continuity ensures that the radial effects 
found here will still be present, and the effects of true lensing will be superimposed. As argued in Q]^], we expect 
caustics to skew the area towards larger values, so that most objects in the universe are demagnified. 

The importance of all this is that it opens up the way for considering the effects of lensing by inhomogeneities on 
large angular-scale number counts and area distances as opposed to limiting discussion to lensing effects on small 
scales. 
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